The ratchet effect and the transporting islands in the chaotic sea 
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We study directed transport in a classical deterministic dissipative system. We consider the 
generic case of mixed phase space and show that large ratchet currents can be generated thanks 
to the presence, in the Hamiltonian limit, of transporting stability islands embedded in the chaotic 
sea. Due to the simultaneous presence of chaos and dissipation the stationary value of the current 
is independent of initial conditions, except for initial states with very small measure. 
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The ratchet effect, that is the possibility of obtaining 
directed transport of particles in the absence of any net 
bias force, is a problem at the heart of statistical me- 
chanics. Ratchet transport in systems at equilibrium is 
forbidden by the second principle of thermodynamics [ij . 
On the other hand, it is possible to overcome this limi- 
tation in systems out of equilibrium, provided all space- 
time symmetries which inhibit directed motion are bro- 
ken [31 • The ratchet phenomenon has recently gained 
renewed interest 0, 0] as a model elucidating the physics 
of molecular motors Moreover, directed transport 
may lead to technological applications at the nanoscale 
Q, including new electron pumps, molecular switches, 
rectifiers and transistors. 

Previous theoretical studies have shown the ratchet ef- 
fect in systems in which noise is absent and its role is 
played by the deterministic chaos induced by the inertial 
term 0]. In particular the origin of current reversal in 
such inertia ratchets has been carefully investigated 0]. 

In spite of these pioneering works, the interrelation 
between the complexity and the rich variety of classi- 
cal chaotic motion in conservative systems and the ap- 
pearance of the ratchet phenomenon when dissipation is 
introduced, is not known. In particular the role of sta- 
ble islands in the mixed phase space structure which is 
generic in nonlinear dynamical systems is not clear Q. 

In this paper, by considering a periodically kicked, dis- 
sipative, inertia ratchet, we show that the generic mixed 
phase space structure of the conservative case may lead 
to a strong ratchet phenomenon when dissipation is in- 
troduced. In particular, for strong dissipation, large cur- 
rents may arise in a short time scale. On the other 
hand, for weak dissipation, large ratchet currents can be 
achieved nearly independently of initial conditions, as a 
result of a beautiful interplay between chaotic diffusion, 
ballistic transport of islands and dissipation. 

The system we study is a particle moving in one di- 
mension [a;e(— oo, cxd)] in a periodic kicked asymmetric 



potential: 



V{x, t) = K[coa{x) + -cos(2a; + (f>)] ^ S{t - mT) , 

(1) 
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where T is the kicking period, which is set to unity in 
this paper. The evolution of the system in one period is 
described by the map 



p = 7P -f K[sin{x) + asin(2a:: - 

X = X + p, 



(2) 



where p is the momentum variable conjugated to x and 
7S[0, 1] is the dissipation parameter, describing a veloc- 
ity proportional damping. The limiting cases 7 = and 
7 = 1 correspond to overdamping and Hamiltonian evo- 
lution, respectively. The spatial symmetry is broken at 
a ^ 0, 7^ riTT, with n integer T2i. However directed 
transport is forbidden in the Hamiltonian limit due to 
time-reversal symmetry [isj . In this limit system ^ ex- 
hibits the typical mixed phase space structure. A central 
point of this paper is to show that even very tiny, almost 
invisible islands, embedded in the chaotic sea play a cru- 
cial role in the generation of large ratchet currents when 
a small dissipation is added. 

In Fig. [T] we plot the ratchet current (p) as a function 
of the dissipation parameter 7 at different times t (the 
discrete time t measures the number of kicks). Here we 
set K = 6.5, a = 0.5, (p — tt/2 and we follow the evolu- 
tion in time of a large number of trajectories whose initial 
conditions are randomly and uniformly distributed in the 
unit cell — tt <p<7r, 0<a;< 27r. Therefore, the ini- 
tial average momentum is (p) = for any 7. Since all 
relevant space-time symmetries Q are broken at 7 7^ 1, 
then a directed current (p) ^ can be generated. The 
dependence of the current on 7 is rather complicated. In 
particular, we can see plateau regions inside which the 
asymptotic ratchet current is independent of 7. We note 
that for strong dissipation (7 < 0.6) the current con- 
verges to its asymptotic value very rapidly and typically 



2 



independently of initial conditions, the dynamics being 
characterized by a single stationary distribution. The 
weakly dissipative regime will be discussed later. 



20 n 10 initial states are 
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FIG. 1: Ensemble average current (p) at difTerent times, for 
a = 0.5. 

In order to understand the behavior of the ratchet cur- 
rent we first perform a linear stability analysis of map ([2]) . 
The fixed points of the map when x is taken modulo 2tt 
are given by 



p* = 21tt, I integer, 

(7 - l)2/7r + K[sm{x*) + acos(2a:*)] = 0, 



(3) 



where for the sake of simplicity we have considered the 
case (j) = tt/2. The stability of these fixed points is de- 
termined by the eigenvalues Ai,2 of the Jacobian 



djp, x) 
dip, x) 



7 K[cos{x) — 2asin(2x)] 
7 1 -I- K[cos{x) — 2asin(2a;)] 



(4) 



The modules |Ai,2| versus x are shown in Fig.[2l^a) for 
different 7. The stable intervals, in which the modules 
of the eigenvalues are both less than 1 , depend only very 
slightly on 7. Therefore we may consider the case 7 = 1, 
for which the stability intervals can be computed analyt- 

1 TT .1 

ically and are given by (arcsin — , — ), (tt — arcsin — , ci), 

4a 2 4a 

(c2, — ) where ci,C2 are the real roots of the equation 

4 

cos(a;) — 2asin(2a;) — . These intervals are shown as 

K 

shadowed regions in Fig.[2l^b). 

As illustrated in Fig. [2jb), the fixed points {x* ,p* — 
2ttI) can be determined graphically. If x* resides in a 
shadowed region, then the fixed point is stable. Clearly 
there are no stable fixed points for ^=0. Instead, for 
any positive integer Z, stable fixed points exist from 7 = 
K{a+-i^) ^ 0.77588 ^ K{l-a) 



2ttI 
0.51725. 
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to 7 = 1 — ■ 
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FIG. 2: (a) Modules of the eigenvalues Xi^2{x) of the sta- 
bility matrix (j4| for different values of 7. (b) Dashed curve: 
/(a::)=A'[sin(2;)+fflcos(2a::)]; solid lines: g{x)={l—^)2lTT for dif- 
ferent I and 7. In the shadowed regions, |Ai,2(2:)| are both less 
than one. 



The bifurcation diagram of Fig. [3] confirms the above 
analytical estimates. The positions of the the main fixed 
point windows {I = 1,2,3) coincide with the plateaus 
(p) — 2tt, 47r, 6tt observed in Fig. [1] At both ends of 
the fixed point windows tangent bifurcations occur which 
correspond to transitions from simple to strange attrac- 
tors. (A similar picture holds for stable periodic orbits, 
even though their positions can not be so easily calcu- 
lated analytically). 
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FIG. 3: Bifurcation diagram: 5 x 10 iterates are drawn after 
a transient of 10^ map steps, starting from 3600 initial states 
drawn from a uniform distribution in the unit cell < x < 2tt, 
—TV < p < n. The analytical expectations for the position of 
fixed point windows at a = 0.5 are also indicated. 

Since the width of the l-th stability window is propor- 
tional to J and ^IZi 7 does not converge, these windows 
must overlap when I is large enough. Therefore, multiple 
attractors (and also periodic orbits) must coexist when 7 
approaches 1 . Their attractive basins cut the phase space 
into many pieces and here one may expect the asymptotic 
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ratchet current to depend, in general, on initial condi- 
tions. The large negative values of (p) that appear at 
long times in Fig. [1] are the consequence of trajectories 
ending up on periodic orbits after a long transient chaotic 
motion. 

For very small dissipation (7 > 0.98), we could only 
observe chaotic motion on a fractal set even though it 
cannot be excluded that this is a transient with a life- 
time much exceeding the total integration time. On the 
other hand this weakly dissipative regime possesses very 
interesting features. Due to chaotic diffusion, the mo- 
mentum probability distribution (Fig. |4]) widens in time 
and eventually, due to dissipation, saturates to a station- 
ary distribution close to a Gaussian. This is clearly seen 
in Fig. [4] in which the most remarkable feature is the 
small peak moving ballistically in the direction of posi- 
tive p. This peak is due to the presence of small stability 
islands in the Hamiltonian limit 7=1 and play a key 
role in the generation of the ratchet current. 




FIG. 4: Snapshots of the momentum probability distribu- 
tion for 1 — 7 = 10~*, starting from 3.6 x 10® initial con- 
ditions randomly and uniformly chosen inside the unit cell 
p £ [—iT,ii),x € [0, 27r). In the inset the two curves drawn 
after long integration times, t = 10^ and t = 5 x 10^, overlap. 
This indicates that the distribution saturates (at least up to 
the integration times explored in our numerical simulations). 

To illustrate this point, we plot in Fig. [5] the evolu- 
tion in time of the first two moments of the probabil- 
ity distribution P{p), and compare two different cases 
corresponding to two different initial conditions: i) uni- 
form distribution inside the unit cell p S [— 7r,7r),a; £ 
[0, 27r), ii) uniform distribution inside the main trans- 
porting stability island. This island, for 7=1, cor- 
responds to ballistic motion in the positive momen- 
tum direction ("accelerator mode"), with the momen- 
tum increased of 47r every 3 map iterations. In- 
deed this island is centered around the periodic orbit 
(xi,pi)«(1.6073, 4.9802) {x2,P2)^ (2.9103,1.303) ^ 
(2^3, P3)«(2. 9103,0) (x4,P4) = (a:i,_pi-l-47r=pi) (here a; 
and p are taken modulo 27r). The area of the island 



w 10~^, that is only 3 x 10~^ of the available phase space. 
Notwithstanding the accelerator mode, the ratchet cur- 
rent avera ges to zero in the Hamiltonian limit: a sum 
rule exists [l05 such that the motion of the islands in the 
direction of positive momentum is balanced by the mo- 
tion of the chaotic sea in the opposite direction. This 
sum rule no longer applies when dissipation is included 
since the latter mixes the sets which are invariant in the 
Hamiltonian limit, namely the islands and the chaotic 
sea. Still, if dissipation is weak, these two sets remain es- 
sentially disconnected for a long time scale tm oc . As 
a consequence, as shown in Fig. [5l for t<tm the ratchet 
current (p)RiO when starting from the entire unit cell, 
while it grows linearly in time when the initial distribu- 
tion is concentrated inside the island. In this latter case, 
as expected, the acceleration of the island is |(27r), that 
is, the island is shifted along p by 2 unit cells every 3 map 
iterations. At t^tm, a momentum high enough is reached 
to allow dissipation to drive particles outside the island. 
Then the motion of these particles becomes chaotic and 
therefore the second moment of the distribution suddenly 
increases, while the first moment decreases due to dissi- 
pation. Correspondingly, for the initial condition inside 
the entire unit cell, the ratchet current starts to increase 
at t~t,„ until the asymptotic value is reached. This pro- 
vides direct numerical evidence that the ratchet current 
is generated thanks to the presence of integrable islands 
in the Hamiltonian limit, the island and the chaotic sea 
being finally mixed by dissipation. 
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FIG. 5: First and second moment of the distribution P{p) as 
a function of time, for two different initial distributions. 

It is interesting to point out that while the ratchet 
current - i.e. the first moment of the distribution P(p) 
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- is determined by a sophisticated dynamical effect, the 
second moment can be understood by means of a simple 
Fokker-Planck equation: 

dP d (\ dP\ d 

where D oc is the diffusion coefficient. The station- 
ary solution to this equation is a Gaussian distribution 



of variance 



D 



which turns out to be in very good 



2(1-7) 

agreement with our numerical data. 

We would like to stress that, thanks to the presence of 
integrable islands in the Hamiltonian limit, large ratchet 
currents can be achieved also for weak dissipation: as 
shown in Fig. [6l we can have {p) > 26. 

A completely different behavior takes place when this 
structure of islands is absent. Consider for example the 
ratchet model: 



generic for systems with spatial and temporal periodicity, 
in that transporting islands are typical for such systems 
in the Hamiltonian limit. 

Finally, we point out that the dissipative ratchet model 
discussed in this paper could be realized by means of cold 
atoms in optical lattices, where it is possible to imple- 
ment the asymmetric pot ential ([1]) [16] and a velocity 
proportional damping |17| . Since kicked systems simi- 
lar to our have been implemented in the deep quantum 
regime [iSj , it appears possible to investigate experimen- 
tally the impact of purely quantum effects such as dy- 
namical localization on the ratchet transport discussed 
in this paper. 

We thank Gabriel Carlo for useful discussions. This 
work is supported in part by an Academic Research 
Grant from MOE and the DSTA under Project Agree- 
ment No. POD0410553. 



p = "fp + K[x — TT — acos(a::)], 
X — X + p. 



(6) 



This map, in the Hamiltonian limit and for a S [0, 1) is 
completely chaotic with no stability islands [3] ■ For this 
map, we have always found as stationary distribution a 
strange attractor, supporting a very weak ratchet cur- 
rent. For instance, at = 1, a = 0.7 we have obtained 
(p) < 0.13 for any value of 7. Basically the ratchet cur- 
rent is only due to the asymmetry of the attracting set 
which, at least for weak dissipation, is also weak. 
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FIG. 6: Stationary ratchet current in the weak dissipation 
regime. The current is averaged over 2 x 10^ map iterations, 
taken after an initial transient of 2 x 10^ map steps. 

In summary, the results presented in this paper show 
that large ratchet currents can be generated in a dissi- 
pative system thanks to the presence, in the Hamilto- 
nian limit, of transporting integrable islands embedded 
in a chaotic sea. This phenomenon leads, due to the 
joint presence of chaos and dissipation, to large ratchet 
currents nearly independently of initial conditions and is 
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